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Abstract. We report new determinations of the radius 
and mass of the pulsating helium-rich, hydrogen-deficient 
star LSS 3184 (BX Cir) using measurements of radial ve- 
locity and angular radius throughout its pulsation cycle. 
Measurements of radial velocity, and thus changes in stel- 
lar radius (Ai?*), were made using Anglo- Australian Tele- 
scope echelle spectra. Hubble Space Telescope ultraviolet 
spectra and ground-based BV photometry were used to 
find temperatures and fluxes throughout the pulsation cy- 
cle. The temperatures and fluxes were used to find the 
angular radius of the star (a). The a, Aa, and AR values 
thus found were used to calculate the mean stellar radius 
(Ri) = 2.31±0.10i? Q . If we use the previously determined 
logg = 3.35 ± 0.1 for LSS 3184 and our radius estimate, 
we find its mass to be M* = 0.42 ± 0.12M©. 
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1. Introduction 

Pulsations in stars provide tools for researchers in sev- 
eral fields of astronomy. They provide standard candles for 
measuring Galactic and extragalactic distances and they 
provide methods for measuring stellar parameters, even 
below the directly observable photosphere. Pulsations in 
hydrogen-deficient, helium-rich stars such as extreme he- 
lium, R Coronae Borealis, and hydrogen-deficient carbon 
stars have not been studied in as much detail as those 
in stars with more 'normal' chemical abundances. While 
this is partly understandable since most pulsating stars 
have normal compositions, studying stars with little or no 
hydrogen is important to allow tests of pulsation theory. 
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Extreme helium stars, as their name would imply, have 
weak or non-existent hydrogen absorption lines and very 
large helium abundances (<; 99 per cent). Two extreme 
helium sta rs, V 652 Her and LSS 3184, are known to pul- 
sate. Saio ( 1993 ) showed that in extreme helium stars with 
temperatures around 2 x 10 5 K, like V652 Her, the re- 
mechanism caused by iron-group (Z-bu mp) o pacity can 
excite the observed pulsations. Saio ( 1994 , 1995 ) predicted 
that LSS 3184 should pulsate because of its l ocatio n in the 
Z-bump instability finger. Kilkenny & Koen ( 1995| ) discov- 
ered that LSS 3184 shows photometric variations with a 
period of about 0.107 d. The similarity of V652 Her and 
LSS 3184 in temperature, surface gravity, and pulsation 
period implies that they are very similar in other physical 
parameters. 



Kilkenny et al. (1999) have recently reported an obser- 
vational analysis of LSS 3184, including a determination of 
its photometric period (0.1066 d). In addition, they used 
medium resolution (A/AA w 4000) spectra to measure 
radial velocity variations and to show that the photomet- 
ric variability is caused by pulsations. Drilling, Jeffery, & 



Heber (1998) reported an analysis of LSS 3184 in which 
they found T cS = 23 300 ± 700 K, l ogg = 3.35 ± 0.1 and 
^h/'t-Hc < 0.00015. Kilkenny et al. (199E) report that us- 



ing the radius they determined and logg from Drilling 
et al. p98| ) gives a mass of 0.15M Q , which is much smaller 
than 0.7M Q , the mass accepted for V652 Her (Lynas-Gray 
et al. 1984), and small enough to imply that some input 
parameter or procedure is in error. 

In this paper we report an analysis of LSS 3184 us- 
ing high resolution optical spectra for radial velocity mea- 
surement and ultraviolet Hubble Space Telescope (HST) 
spectra and ground-based BV photometry for tempera- 
ture, luminosity, and angular radius measurement. The 
new data provide better temperature and angular radius 
estimates and a much cleaner radial velocity curve for the 
star, allowing a more reliable estimate of its radius and 
mass than was possible previously. 



2. Observations and data reduction 
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Fig. 1. Plot showing velocity and phase coverage for LSS 3184 spectroscopic observations, with error bars shown for 
velocity. Note that the velocity here is not corrected for projection effects. 



2.1. AAT visible spectral observations 

Spectra of LSS 3184 were obtained during the nights of 

1996 May 18 and 19 using the University College Lon- 
don Echelle Spectrograph at the 3.9-m Anglo-Australian 
Telescope. Exposure times were between 4 and 5 min- 
utes. Standard iraf packages were used for bias and 
flat field correction, reducing echelle orders to one di- 
mensional spectra, and applying the wavelength scale us- 
ing thorium-argon spectra. As the wavelengths covered 
by adjacent orders overlapped, the spectra covered the 
range 3850-5055 A completely. The spectral resolution 
was A/AA w 48 000. Velocity corrections for Earth's mo- 
tion were found for each exposure using RVCORRECT and 
were applied using DOPCOR. 

2.2. Hubble Space Telescope ultraviolet observations 

Ultraviolet spectra of LSS 3184 were obtained using the 
Faint Object Spectrograph of the HST. Observations were 
made in RAPID mode using the blue detector, the 0'.'86 
square aperture and grating G160L over three orbits on 

1997 February 7. Data files reduced using the most recent 



calibration files were downloaded from the HST archives 
on 1999 October 12. Individual exposures were made ap- 
proximately 19.25 seconds apart. The usable part of the 
spectra covered the range 1150-2500 A. Based on the 
ephemeris of Kilkenny et al. ( |l999| ), observations dur- 
ing the three orbits covered the phase (</>) ranges from 
0.3154-0.6099, 0.9327-0.2376, and 0.5602-0.8651, where 
maximum V magnitude is at <f> = and again at <j) = 1 . 
Thus the star was observed by the HST over 85.5 per cent 
of its pulsational cycle. 



3. Analysis 

3.1. Radial velocity determinations 

The velocity shifts between spectra were measured using 
the cross correlation package FXCOR in iraf. For each 
measurement the velocities found from 28 of the 35 orders 
of the echelle spectra were used to find a weighted average 
velocity. Weights for the averaging were the inverse of the 
velocity errors reported by FXCOR. The unused spectral 
orders either had no strong absorption lines or had strong 
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interstellar lines which did not allow a reliable stellar ve- 
locity determination. 

The procedure took two iterations. In the first iteration 
the best results were obtained by using the sum of all the 
second night's LSS 3184 spectra as the cross correlation 
template for the first night's spectra and vice versa. The 
spectra from both nights were then shifted by the veloci- 
ties thus found and co-added to provide the template for 
the second iteration. The absorption lines in the second 
template were much sharper, as the velocity smearing due 
to the stellar pulsations was effectively removed. A third 
iteration was performed with the spectra shifted by the 
velocities found in the second iteration before co-adding 
to make the template, but the velocities from the third 
iteration were effectively identical to those from the sec- 
ond. The velocities from the second iteration, shifted so 
that the mean velocity is zero, are shown in Fig. |l|. Error 
bars show the weighted standard deviations of the veloc- 
ity from the echelle orders. Gaps are present in the data 
where clouds prevented observations. From the the data in 
Fig. [j], we find that the peak-to-peak radial velocity vari- 
ation is 40.0±0.3 km s _1 , which is larger than 30 km s _1 , 



the value found by Kilkenny et al. (199£). 

Velocity data from the two nights were phased to 
the pulsation cycle using the ephemeris of Kilkenny 
et al. ( |1999[ ): T = 2 449 477.4691(±0.0016), Period = 
0.1065784(±0.0000005)d. The measured radial velocities 
were multiplied by the factor —1.42 to correct for pro- 
jection effects and give the surface velocity through the 
pulsation cycle in the stellar rest frame. The projection 
factor was chosen based on preliminary work by Montahes 
Rodriguez et al . ( |200C| ) (See also Albrow & Cottrell |l994; 
Gautschy 1987; and references therein). We will discuss 
later the effects of choosing a different projection factor. 
A smooth curve (high order polynomial) was fit to the 
velocity data (Fig. ||). The velocity values on this curve 
were used to calculate change in stellar radius (AR+) and 
surface acceleration. The phase bin centers for accelera- 
tion and Ai?* are shifted by half a bin with respect to 
the velocity bins, i.e. the end of an acceleration bin is the 
center of the next velocity bin and the end of a velocity 
bin is the center of the next AR+ bin. 



3.2. Temperature determination 

To find temperature variations through the pulsation cy- 
cle, synthetic spectra were fit to the ground-based B V pho- 



tometry of Kilkenny et al. (1999) and our HST ultraviolet 
spectrophotometry. Before the fitting, the spectra were 
binned in wavelength. No information is lost through the 
binning since the fits are to the shape of the spectrum, not 
to individual lines. The spectra are noisy at short wave- 
lengths, reflecting a drop in detector sensitivity. To avoid 
possible problems caused by this increased noise, only the 
part of the spectra with A > 1270 A was used. 




phase 



Fig. 2. Fit (solid line) to velocity data corrected for pro- 
jection effects (points). The data have been folded over by 
0.2 cycles in Figs. |, |, and §. 



The pulsation period was divided into phase bins and 
the spectra and photometry within each bin were aver- 
aged. No significant difference in the results was found 
with the period divided into 15, 25, 50, or 99 bins. We 
will report the results found for 25 bins. 

A grid of line-blanketed model atmospheres was cal- 
culated under the assumption of plane-parallel geometry, 
hydrostatic equilibrium and local thermodynamic equilib- 
rium using the code Sterne described by Jeffery & Heber 
( |1992| ) and by Drilling et al. Ql998| ). Following the lat- 
ter, we assumed a composition for LSS 3184 given by 
rtH = 0, ?iHo = 0.99 and nc = 0.003, where n represents 
fractional abundance by number, and all other elements 
were assumed to have solar-like relative abundances. The 
grid extended between 15 000 and 30 000 K, log .9 = 3.00 
and 4.00 (cgs units). As will be shown later, changing the 
logg used in the model atmospheres makes only a minor 
difference in the temperatures and angular radii derived. 

In the fitting procedure T c g, angular radius (a), and 
i?B-v were allowed to vary and the downhill simplex pro- 
gram amoeba (Press et 



al. 



1992) was used to find the 



minimum \ difference between the synthetic spectrum 
and the observed spectral and photometric data at each 
phase (Jeffery et al. 2000| ). An example of the fit is shown 
in Fig. |. 

When the extinction -Eb-v was allowed to vary we 
found the mean to be £b-v = 0.239 ± 0.008. Because the 
extinction is not expected to vary with pulsation phase, 
we chose to use £"b-v = 0.24 throughout the cycle and 
did a second iteration allowing only temperature and a to 
vary. 

3.3. Radius determination 

In determining the radius of LSS 3184, we make two as- 
sumptions. First, we assume that the temperature (and 
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Fig. 3. Best fit synthetic spectrum (dashed curve) plotted 
with ultraviolet spectrum (histogram) and BV photome- 
try values. This is for the phase bin centered at 0.580. 



thus a) and radial velocity (and thus Ai?*) were measured 
at approximately the same layer in the stellar atmosphere. 
Second, because we measure a perpendicular to the line 
of sight and Ai?* parallel to the line of sight, we are as- 
suming that the pulsation is spherically symmetric. Thus 
we get Ai?*/i?+ = Aa/a. With these as sump tions we can 
use a modified Baade's method (Baade 1926| ) to find i?*, 
using our previously determined values of a and Ai?* . 

There are two ways we can do this. In the first, we 
choose two phase bins, determine Aa/a and Ai?* between 
them and find i?* = Ai? i a/A«, with a and J?+ defined 
at one of the chosen phase bins. This method has the 
disadvantage of using data from only two phase bins and 
thus ignoring additional information available from the 
rest of the pulsation cycle. In the second method we plot 
Aa/a versus AR+. The slope of a linear fit to the data 
points is then We use the second method. 



4. Results and discussion 

The velocities, effective temperatures, and angular radii 
measured through the pulsation cycle of LSS 3184 are 
listed in Table 1. The velocities are not corrected for pro- 
jection effects and are reported for the center of the phase 
bins. Temperatures and angular radii are reported for the 
beginning of the phase bins. Data are missing where HST 
observations did not fall in the affected phase bins. Typical 
uncertainties are indicated. 

The surface acceleration, surface velocity, and change 
in radius determined for LSS 3184 using the AAT spectra 
are shown in Fig. [|. Note that the surface velocity reported 
here is the velocity measured from the spectra multiplied 
by —1.42 to correct for projection effects and make posi- 
tive velocity be away from the star's center. The change 
in radius has been set so that it is zero at photometric 
maximum, = 0. 



Table 1. Velocity, temperature, and angu- 
lar radius through the pulsation cycle of 
LSS 3184. Velocity is not corrected for pro- 
jection effects. 
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Fig. H shows the integrated flux between 1270 and 
2508 A and the temperature and angular radius deter- 
mined by fitting synthetic spectra to HST UV spectra 
and ground-based B V photometry through the pulsation 
cycle. 

Aa/ao is plotted against AR* in Fig. ^. The '0' sub- 
script indicates that the reference bin is at <f> = 0. As 
is seen in Figures ^ and [sj, the shapes of the a and Ai?* 
curves are not identical, which is the cause of the non- 
linear, looped shape of the Aa/a versus AR± curve. If 
the a and Ai?* curves are normalized and placed on the 
same plot (Fig. |?]), the differences in the c urves are more 
noticeable. Figure 12 from Kilkenny et al. ( |1999| ) shows a 
similar loop on the right hand side of the angular radius 
versus stellar linear radius plot. Our work cannot be used 
as an independent confirmation of the non-linear angular 
versus stellar radius curve, however, as we have used the 
same BV photometry as Kilkenny et al. 

The non-linear shape makes determining the radius 
more difficult. It raises questions about the assumption 
made in using Baade's method to find R+, that AR+ 
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Fig. 4. Surface acceleration, surface velocity, and change 
in radius of LSS 3184 through its pulsation cycle. 



and Aa are measurements of the same quantity, with 
Aa decreased by a factor proportional to the distance to 
LSS 3184. One possible explanation of the discrepancy is 
that we are measuring Ai?* and a at different layers in 
the atmosphere. Ai?* is measured using optical spectra, 
while a is measured using ultraviolet spectra and optical 
photometry. It is possible that the layer where the optical 
lines are formed expands and contracts a bit differently 
than the layer where the ultraviolet continuum is formed. 
It is also possible that the star and/or its pulsations are 
not spherically symmetric, as might occur if the star were 
flattened by rotation, so that the measured Aa, which 
is perpendicular to the line of sight, and Ai?*, which is 
parallel to the line of sight, act differently. Further, the 
atmosphere is not static. The atmosphere's temperature 
is constantly changing, and the pulsating layers undergo 
a substantial compression at minimum radius, as shown 
by the spike in the acceleration in Fig ^, which may cause 
nonadiabatic effects. So the temperature may be acting 
differently in the expanding part of the phase than in the 
contracting part. 



Fig. 5. Flux (Fi27o-2508)j effective temperature, and an- 
gular radius of LSS 3184 through its pulsation cycle. The 
curves simply connect the data points and are not fits to 
the data. 



However, it is encouraging that the slopes of the ex- 
panding (the upper half of the loop in Fig. ^) and con- 
tracting (lower half of the loop) parts of the Aa/a versus 
AR± curve are not too different. Aa and AR* are still 
correlated. 

If we find the slope of the curve using a least squares 
fit to all of the data points (dashed line in Fig. ^J) then 
the radius at photometric maximum (0 = 0), the inverse 
of the slope, is 2.27 ± 0.10i?o, where the uncertainty is 
derived from the standard error in the least squares fit 
to the slope. The average Ai?* over the pulsation cycle, 
with AR* = at = 0, is 28 190 km, or 0.04i? o . So if 
we use all the Aa/ao versus Ai?* data points, we find the 
mean R* to be (i?*) = 2.31i?Q. This is larger than the 
1.35i?0 mean radius found by Kilkenny et al. (1999) for 



LSS 3184 and is closer to the 1.91i?n_mean radius found 



for V652 Her by Lynas-Gray et al. (|1984|). 

We tested the effects of using only a portion of the 
Aa/a versus Ai?* curve to determine i?*, though we have 
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This is larger than 0.15M©, the mass found by Kilkenny 
et al. Q1999D for LSS 3184 and closer to 0.7i°t, the esti- 



Fig. 6. Aa/ao versus Aft* through the pulsation cycle of 
LSS 3184. Points are connected in order of phase with 
<fr = at the origin. The dashed line is a linear least 
squares fit to all data points. Error bars for Aft* are 
smaller than the symbols. Numbers next to symbols in- 
dicate the corresponding pulsation phase. 




Fig. 7 

sation 



, a (solid curve) and Aft (dashed curve) versus pul- 
phase. 



no reason to reject any particular data points. If we use 
the points on the upper part of the curve, with phase 
0.16 < 4> < 0.56, then we find (ft*) = 2.16 ± 0.07ft©. If 
we use the points on the lower part of the curve, ignoring 
the points to the left of the change in slope at Aft* « 
10 000 km, so that 0.60 < <f> < 0.96, we find (ft*) = 1.75 ± 
0.08ft©. If we use only the points to the right of Aft* = 
km then we find (ft*) = 1.97 ± 0.10ft©. In all cases, 
even the extreme one where we reject all points but those 
on the lower part of the curve, our (ft*) is larger than 
ft* = 1.35 ± 0.15ft©, the value found by Kilkenny et al. 



(|1999D 



mated mass of V652 Her (Lynas-Gray et al. 1984). Our 
larger estimate results mainly from the larger peak-to- 
peak range of radial velocities we measured. It is likely 
that the smaller velocity amplitude measured by Kilkenny 
et al. resulted from a combination of the lower spectral res- 
olution and the lower signal to noise in the spectra they 
used. This meant that even the best cross correlation tem- 
plate still contained some velocity broadening, thus dilut- 
ing the velocity amplitude. 

There are several sources of uncertainty in our mass 
determination. Uncertainty in the logg used is a major 
contributor. We note that it is impossible to determine 
the pulsation phase of LSS 3184 when the spectrum used 
in the analysis of Drilling et al. ( 199S| ) was taken in 1985, 
as P is unknown. They estimated the temperature of 
LSS 3184 at 23 300±700 K. Our data yield a similar value, 
(T cff ) = 23 230 K, if we use £ B -v = 0.27, the extinction 
they used. However, because temperature, line strengths, 
and other parameters presumably varied throughout the 
1-hour exposure for their spectrum, it is possible that the 
errors in their temperature and gravity determinations are 
larger than the formal errors quoted. If we assume that 
LSS 3184 has the same logg as V652 Her, logg = 3.68 
(Jeffery et al. 1999), instead of logg = 3.35, then we find 
that it has M* = 0.92M Q 

As mentioned earlier, the model atmospheres used 
to calculate temperature and angular radius from HST 
UV spectra and ground-based BV photometry assumed 
logg = 3.50. The assumed gravity enters into the stellar 
parameter calculations in the stellar atmospheres used and 
in calculating the mass from the radius. Changing logg in 
the model atmospheres by 0.50 dex (to 3.00 or 4.00), but 
using logg = 3.35 to calculate the mass as before, changes 
the (T c g), (a), and (ft*) derived by about 1 per cent, and 
thus has only a small effect (~ 3 per cent) on the mass 
derived. The mass finally derived is proportional to g. 

Uncertainty in the extinction also adds uncertainty to 
the stellar parameters calculated. Using £b-v = 0.25 in- 
stead of 0.24 increases the derived (T G g) by 510 K, or 
about 2 per cent, but changes (a) by only 0.2 per cent, 
and thus has a very minor effect on the derived stellar 
radius and mass. 

The derived radius varies proportionally with the pro- 
jection factor used to transform the measured stellar ra- 
dial velocities into surface velocities. For example, using a 



Drilling et al. (1998) estimated logg = 3.35 ± 0.1 for 
LSS 3184. If we use this with our (ft*) = 2.31ft© in the 
formula g = GMR~ 2 , we find M* = 0.42 ± 0.12M Q . 



projection factor of 1.31, as Lynas-Gray et al. (1984) used 
in their analysis of V652 Her, instead of 1.42 would give 
(ft*) = 2.12ft Q instead of (ft*) = 2.31ft©. The smaller ra- 
dius would give M* = 0.37M© instead of M* = 0.42M©. 

5. Conclusions 

In this paper we report new determinations of the radius 
((ft*) = 2.31±0.10ft©) and mass (M* = 0.42±0.12M©) of 
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the pulsating hydrogen-deficient star LSS 3184. The radial 
velocity data and temperature measurements are more re- 
liable than those used for previous radius and mass deter- 
minations for this star, so our estimates are likely to be 
closer to the star's actual parameters. Further improve- 
ments can be accomplished by improving the log g and 
chemical composition estimates for LSS 3184. In addition, 
there needs to be further study of the difference between 
the shapes of the angular radius and physical radius curves 
based on temperature and flux measurements, and radial 
velocity measurements, respectively. It is important to de- 
termine if the difference in shapes is a sign that some of 
our data or a step in our analysis is flawed, or if there is 
something happening in the star itself that causes the two 
measures of radius to behave differently with pulsational 
phase. 
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